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NOTATION 


B 

B 


CB 

CG 

K,  M,  N 

kb,  mb,  nb 
Kg,  Mq,  hG 


Magnitude  of  B 

Resultant  force  caused  by  static  pressure  of  fluid  acting  on 
boundary  of  body 

Center  of  bouyancy  (the  point  where  B  appears  to  act) 

Center  of  mass  (the  point  where  W  appears  to  act) 

Components  of  resultant  moment  of  B  and  IV  about  0  in  the 
directions  of  the  x,  y,  2  axes  4 

Moments  of  B  about  the  x,  y,  2  axes 
Moments  of  W  about  the  x,  y,  2  axes 


0 

Q 

Qb 

Qg 

rR 


Origin  for  the  x,  y,  2  axes 
Combined  moment  of  B  and  W  about  0 
Moment  of  B  about  0 

Moment  of  W  about  0 

Position  vector  of  CB  with  respect  to  0 


r 


G 


Position  vector  of  CG  with  respect  to  0 


IV 

w 

X,  Y,  Z, 

Xq>  etc. 


Magnitude  of  W 

Resultant  force  caused  by  gravitational  attraction  on  mass  of 
body.  Water  inside  the  boundary  of  the  body,  whether  in  hard 
tanks  or  free-flooding,  normally  is  considered  part  of  the  mass. 

Components  of  the  resultant  of  B  and  W  in  the  directions  of 
the  x,  y,  2  axes 

6X  dX 

Partial  derivatives  —  ,  —  ,  etc. 

d6  dcf> 

Initial  value  of  X  force 


r,  y,  z 


Coordinate  axes  fixed  in  the  body  (See  DTMB  Report  1319 x) 


xn’  Vr'  zr 


Coordinates  of  rB  in  the  x,  y,  2  axes 
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References  are 


!;ited  on  page  9  . 


Ill 


XC>  Vc*  2g 


Coordinates  of  rc  in  the  x,  y,  z  axes 


y,,  a. 


A 


ip,  d ,  <t> 


!^o’  ^o> 


A  set  of  orthogonal,  right-hand,  Cartesian  coordinate  axes 
fixed  in  space  (The  z j  axis  is  assumed  fixed  in  the  direction 
of  gravitational  attraction.) 

Prefix  indicating  an  incremental  change  in  a  quantity 

Euler  angles  giving  the  orientation  of  the  x,  y,  z  axes  with 
respect  to  the  xJ}  y j,  z{  axes  (See  Reference  1) 

Initial  values  of  0,  6 ,  <f> 
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ABSTRACT 


Complete  expressions  are  given  for  the  forces  and  moments  acting  on  a 
submerged  body  because  of  gravitational  attraction  on  the  mass  of  the  body  and 
buoyancy  force  acting  on  the  body  boundary.  The  variations  of  the  forces  and 
moments  with  changes  in  angular  attitude  of  the  body  are  discussed  in  detail. 
The  relationship  of  these  variations  to  the  phenomenon  of  metacentric  stability 
is  described. 


INTRODUCTION 

Frequently  the  motion  of  an  underwater  body  is  treated  by  considering  that  the  body  is 
ballasted  until  it  is  in  neutral  buoyancy,  that  is,  the  mass  of  the  body  is  adjusted  to  equal  the 
displacement  of  the  body.  Coupled  with  such  an  approach,  it  is  usually  assumed  that  the  loca¬ 
tion  of  the  center  of  mass  simultaneously  is  adjusted  to  a  position  directly  under  the  center  of 
buoyancy.  Such  a  body  is  said  to  be  in  neutral  buoyancy  and  level  trim. 

No  great  complication  is  introduced  in  the  equations  of  motion  if  one  assumes  that  the 
body  is  neither  in  level  trim  nor  in  neutral  buoyancy-in  fact,  such  an  assumption  is  usually 
more  representative  of  the  true  condition,  an  out-of-trim  state  either  by  design  or  because  of 
leakage.  As  an  interesting  aside,  leakage  may  be  of  air  at  a  ballast  tank  blow  valve,  thus 
causing  the  body  to  be  “light.” 

The  treatments  of  the  out-of-trim  motion  of  an  underwater  body  which  exist  in  the  liter¬ 
ature  usually  employ  expressions  which  have  been  “Iinerarized”  to  some  degree.  The  author 
felt  that  complete,  nonlinear,  expressions  for  the  gravitational  and  buoyancy  terms  might  be 
useful  because  such  expressions  would  be  exact  regardless  of  the  magnitude  of  the  angular 
departure  from  “even  keel.”  Furthermore  the  proposed  expressions  would  involve  no  restric¬ 
tions  as  to  the  locations  of  the  center  of  buoyancy  and  center  of  mass  with  respect  to  the  ori- 
gin  of  the  body  axes.  Such  exact  expressions  can  readily  be  linearized  to  the  extent  warranted 
by  a  particular  problem.  The  simplification  of  the  basic  expressions  in  the  case  of  a  plane  of 
symmetry  is  illustrated. 

The  material  reported  here  is  part  of  a  broader  study  of  the  general  equations  of  motion 
of  a  body  in  a  fluid  being  carried  out  under  the  General  Hydromechanics  Research  Program  at 
the  David  Taylor  Model  Basin.  Other  related  work  is  contained  in  References  2  and  3. 

GENERAL  CONSIDERATIONS 

The  phenomenon  described  in  the  field  of  naval  architecture  either  as  static  stability 
or,  preferably,  as  metacentric  stability  results  from  the  interaction  of  two  external  forces  act¬ 
ing  on  a  body  that  has  part  or  all  of  its  outer  envelope  in  contact  with  a  fluid.  These  forces 
are  called  the  gravitational  force  and  buoyancy  force.  The  importance  of  metacentric  stability 


in  the  dynamics  of  underwater  bodies  has  been  recognized  for  some  time.4  The  equally  im-  ’ 
portant  metacentric  stability  of  a  surface  ship  is  a  slightly  more  involved  concept  and  lies 
outside  the  scope  of  this  paper. 

BUOYANCY  FORCE 

The  buoyancy  force  frequently  is  called  displacement.  It  is  a  function  of  only  the  form 
of  the  outer  envelope  of  the  body  and  the  presence  of  a  static  pressure  gradient  in  the  fluid. 

The  static  pressure  gradient  exists  only  because  of  the  gravitational  attraction  on  the  mass 
elements  of  the  fluid. 

Surfaces  of  equai  static  pressure  may  be  visualized  as  horizontal  laminations  in  the 
fluid,  with  the  pressure  increasing  as  one  goes  deeper  because  each  lamination  is  supporting 
the  weight  of  all  the  layers  of  fluid  above  it.  Thus  the  direction  of  maximum  static  pressure 
gradient  is  vertically  down. 

Any  body  immersed  in  the  fluid  will  have  some  portion  that  penetrates  into  the  deeper 
laminations  and  hence  will  experience  a  higher  static  pressure  on  this  lower  portion.  It  fol¬ 
lows  that  the  body  is  acted  upon  by  a  net  buoyancy  force  upward,  with  a  line  of  action  directly 
opposite  to  the  direction  of  gravitational  attraction. 

if  the  fluid  and  the  body  are  incompressible,  the  magnitude  of  the  buoyancy  force  is 
directly  proportional  to  the  volume  of  fluid  displaced  by  the  body  and  is  independent  of  the 
angular  attitude  or  depth  of  the  body.  Conceptually,  the  buoyancy  force  may  then  be  con¬ 
sidered  to  act  at  a  fixed  point  whose  location  is  called  the  center  of  buoyancy.  The  position 
of  the  center  of  buoyancy  is  governed  solely  by  the  shape  of  the  body-fluid  boundary. 

Generally  speaking,  hull  compressibility  will  not  cause  significant  migrations  of  the 
center  of  buoyancy  of  an  underwater  body  but  may  account  for  sensible  changes  in  the  magni¬ 
tude  of  the  buoyancy  force.  The  compressibility  of  w$ter  is  considered  of  no  practical  im¬ 
portance  unless  very  great  depth  changes  are  envisioned;  density  is  of  the  order  of  5  percent 
higher  in  the  greatest  ocean  depths  (35,000  feet)  than  on  the  surface.  The  “density  layers” 
caused  by  abrupt  changes  in  salinity  or  temperature  normally  are  much  more  significant. 

GRAVITATIONAL  FORCE 

Gravitational  attraction  acts  on  all  the  mass  elements  of  the  body.  The  distribution  of 
mass  within  the  body  determines  the  location  of  a  point  called  the  center  of  mass,  about  which 
the  body  will  balance  in  a  uniform  gravitational  field.  Strictly  speaking,  the  integrated  effect 
of  gravitational  attraction  on  all  the  mass  elements  of  the  body  gives  rise  to  a  net  gravitational 
force  which  may  be  assumed  to  act  at  a  point  called  the  center  of  gravity.  The  center  of  gravity 
will  differ  from  the  center  of  mass  only  if  the  gravitational  field  is  nonuniform. 
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In  relation  to  the  probable  size  of  underwater  bodies,  the  earth’s  gravitational  field 
may  be  considered  uniform;  furthermore,  the  center  of  mass  is  of  special  importance  if  accel¬ 
erated  motions  of  the  body  occur,  hence  the  preference  for  that  term,  vis-a-vis  center  of  grav¬ 
ity,  in  discussions  of  dynamics. 

In  typical  underwater  bodies,  the  magnitude  of  the  gravitational  force  and  the  location 

of  its  point  of  action  are  adjustable  by  ballasting;  i.e.,  the  mass  and  center  of  mass  are 
variable. 

MET  ACENTRIC  STABILITY 

The  existence  of  metacentric  stability  depends  upon  the  failure  of  the  center  of  mass 
and  the  center  of  buoyancy  to  coincide.  Assume  a  body  submerged  in  water,  motionless,  and 
on  “even  keel.”  Furthermore,  assume  that  the  body  has  been  ballasted  so  that  the  center  of 
mass  is  directly  below  the  center  of  buoyancy  and  that  the  mass  has  been  adjusted  to  produce 
a  gravity  force  (hereafter  called  weight)  just  equal  to  the  buoyancy  force.  Such  a  body  will 
remain  motionless  at  its  present  depth.  If  the  body  then  is  disturbed  in  either  pitch  or  roll 
(and  the  disturbing  cause  subsequently  removed),  the  body  will  return  to  its  initial  attitude. 
Such  a  body  is  said  to  possess  the  characteristic  of  metacentric  stability.  The  stability 
exists  because  any  angular  displacement  about  any  axis  (except  the  vertical  line  containing 
the  center  of  mass  and  center  of  buoyancy)  shifts  the  lines  of  action  of  the  two  forces  to  two 
parallel  locations  such  that  a  couple  is  produced  tending  to  decrease  the  angular  displacement. 

Two  important  variations  on  the  preceding  situation  can  be  pictured  In  the  first,  the 
body  is  in  neutral  buoyancy  but  not  in  level  trim.  As  a  consequence,  the  center  of  buoyancy 
is  not  above  the  center  of  mass.  When  the  body  is  released  in  the  stated  initial  conditions, 
it  will  rotate  to  a  new  attitude  where  the  two  centers  do  lie  on  a  common  vertical  line;  in  the 
new  attitude,  the  configuration  will  have  metacentric  stability,  as  described  before.  Note 
that  a  net  initial  moment  exists  which  tends  to  rotate  the  body  to  its  new  attitude. 

In  the  second  important  variation,  the  center  of  buoyancy  is  above  the  center  of  mass 
but  the  weight  and  buoyancy  are  not  equal.  There  is  a  net  initial  vertical  force  present, there¬ 
fore,  which  will  accelerate  the  body  when  it  is  released  under  the  stipulated  initial  conditions. 
For  discussion,  assume  that  the  weight  is  greater  than  the  buoyancy.  Then  a  fraction  of  the 
weight  equal  to  the  buoyancy  can  be  paired  with  the  buoyancy  to  create  the  conditions  for 
metacentric  stability  described  in  the  first  paragraph  of  this  section.  The  remainder  of  the 
weight  constitutes  an  accelerating  force  at  the  center  of  mass  which  can  exist  simultaneously 
with  the  conditions  for  metacentric  stability. 

A  third  variation  that  combines  the  cases  described  in  the  two  preceding  paragraphs 
does  not  introduce  any  new  concepts. 
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SUMMARIZING  REMARKS 

The  preceding  discussion  has  sought  to  develop  two  principal  points  in  the  general 
case  of  a  submerged  body  not  in  neutral  buoyancy,  namely:  (1)  there  is  an  initial  net  vertical 
force  and  a  net  moment,  at  a  given  initial  attitude,  and  (2)  with  changes  in  angular  attitude 
there  are  changes  in  the  moment,  represented  by  the  phenomenon  of  metacentrie  stability. 
Because  a  net  vertical  force  exists,  its  presence  leads  to  force  changes  along  the  body  axes 
if  the  angular  attitude  changes. 

DERIVATION  OF  GENERAL  EXPRESSIONS 

In  this  section,  the  complete  expressions  will  be  given  first  for  the  force  and  moment 
on  an  out-of-trim  submerged  body.  A  presentation  of  what  are  called  “attitude  derivatives” 
will  then  be  given.  These  derivatives  describe  the  rate-of-change  of  force  and  moment  with 
change  in  angular  attitude. 

EXPRESSIONS  FOR  FORCE  AND  MOMENT 

Assume  a  set  of  orthogonal,  right-hand,  Cartesian  coordinate  axes  xJt  r/j,  z,  fixed  in 
space,  with  the  axis  aligned  in  the  direction  of  gravitational  attraction.  This  set  will  bo 
called  inertial  axes.  Also  assume  a  submerged  body  with  a  similar  set  of  axes  x,  y,  z  fixed 
in  the  body  at  the  point  O.  Let  0,  9,  0  be  Euler  angles  fixing  the  angular  orientation  x,  y,  z 
with  respect  to  xJt  yJf  zf  (see  Reference  1). 

Let  B  be  a  vector  representing  the  buoyancy  force  on  the  body,  and  W  a  vector  repre¬ 
senting  the  force  of  gravitational  attraction  on  the  mass  of  the  body.  Magnitudes  of  the  vec¬ 
tors  are  represented  by  omitting  the  bars  over  the  letters.  Then 

W  m  mg  [l] 

where  m  is  the  mass  of  the  body  and  g  is  the  acceleration  of  gravity. 

Consider  now  the  resultant  of  the  buoyancy  and  weight  forces,  the  vector  sum  of  W  and 
B.  The  magnitude  of  the  resultant  is  W  -B  because  B  is  opposite  in  direction  to  W.  The 
selection  for  the  orientation  of  the  inertial  axes  was  such  that  the  resultant  has  a  component 
only  along  the  zf  axis,  and  its  amount  is  W  -B.  The  resultant  also  has  components  along 
the  x,  y,  z  axes,  which  are  represented  by  X,  Y ,  Z,  respectively. 

The  following  table  of  direction  cosines  can  be  used  to  obtain  the  omponents  X,  Y,  Z 
in  terms  of  the  components  that  the  resultant  has  in  the  inertial  axes. 
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W  -  B 


cos  0  cos  0 


sin  0  cos  0 


-sin  0 


-sin  0  cos  0  +  cos  0  sin  0  sin  0 


cos  0  cos  0  +  sin  0  sin  0  sin  0 


cos  0  sin  0 


cos  0  sin  0  cos  0  +  sin  0  sin  0 


sin  0  sin  0  cos  0  -  ccs  0  sin  0 


cos  0  cos  0 


Thus 


A  «  -  (iy  -  B)  sin  0 
y  =  (W  -  6)  cos  0  sin  0 
Z  ■  (W  -  B)  cos  0  cos  0 

Assume  that  5  has  a  point  of  action  at  the  center  of  buoyancy  CB.  L.9t"7^  be  a  vector 
giving  the  position  of  CB  with  respect  to  0.  Similarly,  let  W  act  at  the  center  of  mass  CG , 
and  let  rG  be  its  position  vector.  Then  B  has  a  moment  about  0 ,  represented  by  the  symbol 
Qg ,  and  given  by  the  vector  cross  product 

Qf)  m  rB  *  ^  [3] 

Let  Qb  have  components  Kg,  Mg,  N B  about  the  x,  y,  a  axes,  respectively.  Then  from 
Equation  [3] 


KB  m  (-  B  cos  0  cos  0)  yB  <■  ( B  cos  0  sin  0)  ag 
Mg  m  ( B  sin  0)  Zg  +  (B  cos  0  cos  0)  x B 
Ng  m  (-5  cos  0  sin  0)xg  -  ( B  sin  0)yg 

where  xB,  yB,  zg  are  components  of  in  the  directions  of  the  x,  y,  z  axes. 

The  weight  force  vector  W  also  has  a  moment  about  0  according  to 

Qg  "  rG  x  ^  [5] 

where  Qc  is  the  moment  about  0. 

Let  rG  have  components  xG>  yG,  zG  in  the  directions  of  the  x,  y,  z  axes,  and  take 
K G,  M G,  Ng as  components  of  Qc  about  those  axes.  Then  Equation  [5]  can  be  separated  into 

Kc  -  (iy  cos  0  cos  0)  yG  -  (If  cos  0  sin  <f>)  zG 
Mg  -  (-iy  sin  0)  a G  -  (iy  cos  0  co 3  $)xG 
f^G  =  (iy  cos  0  sin  <f>)xG  +  (iy  sin  0)yG 
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Let  Q  be  a  vector  representing  the  moment  about  0  of  the  combined  effects  of  B  and  W , 
and  let  K,  M,  N  be  components  of  Q  about  the  x,  y,  z  axes.  Then  from  Equations  [4]  and  [6] 

K  ■  cos  8  [(yclV  -yB B)  cos  <f>  -  (zGW  -sbB)  sin  <f>  ] 

M  =  -  (xqW -xbB)  cos  8  cos  <f>  -  (zGW  -z  bB)  sin  8 

N  a  (xqW -xbB)  cos  8  sin  <f>  +  {VqW  -  yB&)  sin  0 

ATTITUDE  DERIVATIVES 

Assume  that  a  submerged  body  is  at  some  initial  angular  attitude  specified  by  ^Q,  8Q, 
0O.  Then  differentiation  of  Equations  [2]  and  [7]  with  respect  to  8  or  <f>  will  supply  the  fol¬ 
lowing  series  of  complete  expressions  for  the  derivatives.  (Note  that  Equations  [2]  and  [7] 

are  not  functions  of  \p.) 


Xq  =  -  (W  -B)  cos  dQ 

V° 

Yq  a  -  (W  -B)  sin  dQ  sin  0O 
-(*-*)  cos  8q  cos  cf>Q 
Zq  «  -  (W  -B)  sin  dQ  cos  <f>0 
Z<f>B  -(W-B)  cos  dQ  sin  <£0 

k8  "  ~ sin  6o  [(yc^-ysfi) cos  ■( zgw~zbb ) sin  *o] 

K<t>  -  - cos  e0  [  (zcw-zbb)  cos  +(ycw-yBBy> sin  0o] 

Mq  -  (2 GW  -  z bB)  cos  +(xqW  -xbB)  sin  dQ  cos  <60 

%  *  ( xGW~xBB >  cos  d0  sin  ^0 

Ne  =  ( VGW~yDB )  cos  90  sin  °0  sin  ^0 

^  ■  ( xgw~xbb )  cos  0o  cos  ^0 


[81 


Where  Xq  = 


aeT 


etc. 


EFFECTS  OF  LARGE  PERTURBATIONS  ON  ATTITUDE  DERIVATIVES 

Note  that  Equations  [8]  are  partial  derivatives  and,  although  they  are  valid  for  initial 
values  \pQ,  dQ,  <£q  however  large,  they  apply  only  to  snail  incremental  changes  A 0,  A<£  from 
those  initial  conditions. 
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If  the  perturbations  from  the  initial  conditions  are  large,  such  that  the  approximations 


si;>,  Ad  »  0 
cos  Ad  a  1 
sin  A<£  =  0 

COS  A0  a  1 


>  [9] 


J 


are  not  valid,  Equations  [8]  must  be  replaced  by  more  cumbersome  expressions.  For  example, 
let  A0  be  the  initial  value  of  A  force  due  to  an  initial  value  of  dQ  so 

XQ  =  -(W-B)  sin  dQ  [10] 

from  Equations  [2].  After  a  large  perturbation  Ad,  there  will  be  an  attendant  change  in  X  force 
of  amount  AA,  whereby 


AQ  +  AA  =»  -  (W  -B)  sin  (dQ  + Ad)  =  -  (W-B)  (sin  dQ  cos  Ad  +  cos  d0  sin  Ad)  [11] 

Subtracting  Equation  [10]  from  Equation  [11],  the  change  in  A  force  due  to  a  change  in  d  is 

AA  =  -  (W  -  B)  cos  dQ  sin  Ad  +  (IP  -  B)  sin  dQ  (1  -  cos  Ad)  [12] 

Note  the  comparison  of  Equation  [12]  with  the  first  member  of  Equations  [8], 

A  similar  treatment  can  be  given  the  remainder  of  Equations  [2]  and  [7].  As  an 
illustration 


AF  *  -(iy-B)  sin  dQ  sin  <f>Q  sin  Ad  cos  A<£ 

+  (!K-fi)  cos  0Q  cos  <£0  cos  Ad  sin  A<£ 

-(&’-£?)  sin  dQ  cos  0O  sin  Ad  sin  A <f> 

-(W-B)  cos  dQ  sin  <f>Q  (1-cos  Ad  cos  A<£)  [13] 


which  may  be  compared  with  the  third  and  fourth  members  of  Equations  [8].  Equations  of  the 
form  of  [12]  and  [13]  hold  for  perturbations  however  large. 

The  perturbations  may  be  moderately  large  where  the  approximations  of  Equations  [9] 
do  not  hold  but  the  assumptions  can  be  made 

sin  Ad  »  Ad 
cos  Ad  ■  1 
sin  A<£  a  A<£ 
cos  A<£  a  1 
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These  approximations  are  valid  to  within  1  percent  for  angles  up  to  8.1  degrees.  Then  Equa- 
tion  [12]  can  be  replaced  by 


£aX  »  -(W-B)  (cos  $0)  A0 


[15] 


and  Equation  [13]  by 

AF  «  -(W-B)  j^(sin  dQ  sin  <^)A0-(cos  0Q  cos  <^0)A^  +  (sin  dQ  cos  0^)  AtfA^j  [16] 

Note  that  A0A0  is  less  than  0.02,  in  radian  measure,  when  the  two  perturbations  are  each 
8  degrees.  The  term  involving  A0A 0  in  Equation  [16],  therefore,  is  of  the  same  order  as  the 
term  involving  1-cos  A0  cos  A 0,  which  was  dropped  from  Equation  [13]  as  a  result  of  applying 
Equations  [14].  If  the  term  involving  A0A0  is  also  dropped,  Equation  [16]  becomes 

AF  -  -(1F-B)  £(sin  sin  0o)A0-(cos  0Q  cos  0^)  A0^  [17] 

Comparison  will  show  a  mathematical  equivalence  of  Equation  [15]  with  the  first  two 
members,  and  of  Equation  [17]  with  the  third  and  fourth  members,  of  Equations  [8].  The  con¬ 
clusion  is  made,  therefore  that  Equations  [8]  are  suitable  for  perturbations  of  up  to  8  degrees 
in  6  and  0  from  any  arbitrary  initial  conditions  of  these  quantities,  provided  that  errors  of 
less  than  2  percent  are  acceptable.  For  larger  excursions,  or  a  more  stringent  accuracy  re¬ 
quirement,  expressions  of  the  form  of  Equations  [12]  and  [13]  should  be  used. 

EFFECT  OF  PLANE  OF  SYMMETRY  ON  ATTITUDE  DERIVATIVES 

If  the  body  possesses  a  common  plane  of  symmetry  for  both  its  geometric  shape  and  its 
mass  distribution,  the  attitude  derivatives  can  be  given  in  simpler  form.  Assume  that  the 
y  axis  is  normal  to  the  plane  of  symmetry  and  that  the  origin  for  the  x,  y,  z  axes  lies  in  the 
plane  of  symmetry.  Then 

y  b  m  y  g  ~  ®  usi 

and  Equations  [8]  can  be  replaced  by 


Xq  =  -(W-B)  cos  6q 
*0-0 

y q  -  -(W-B)  sin  0O  sin  0O 
y0  -  (W~B)  cos  eQ  cos  0O 


[19] 


* 


Zq  a  -(If'-B)  sin  'JQ  cos  $0 
Z ^  ■  -( W-B )  cos  0C  sin  <pQ 
K q  =  (zqW -ZgB)  sin  dQ  sin  <f>Q 

K<f>  3  "( zgW~zBB )  cos  9o  cos 

A/^  a  -(2 GW -z bB)  cos  60+{xqW-xbB)  sin  dQ  cos  <£0 

^  =  {xgW~xBB)  cos  0o  sin 
Sq  *  -(ar^W'-jr^fi)  sin  0Q  sin  <f>Q 

N4>  m  ( xcW~xBB )  cos  0o  cos  $0 
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